Scalar perturbations of higher dimensional rotating and ultra-spinning black holes 
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We investigate the stability of higher dimensional rotating black holes against scalar perturbations. 
In particular, we make a thorough numerical and analytical analysis of six-dimensional black holes, 
not only in the low rotation regime but in the high rotation regime as well. Our results suggest 
that higher dimensional Kerr black holes are stable against scalar perturbations, even in the ultra- 
spinning regime. 



PACS numbers: 



I. INTRODUCTION 



Exact solutions to Einstein equations are extremely 
useful, specially if they describe simple yet physically 
attainable systems. Indeed, take for example the fa- 
mous Schwarzschild metric: with this exact solution 
at hand, describing the geometry outside a spherically 
symmetric distribution of matter, one was able to com- 
pute the deflection of light as it passes near the Sun 
(and to match the theoretical prediction against the ob- 
servational data), thereby giving a strong support to 
Einstein's theory. We now know that the outside ge- 
ometry of many astrophysical objects is well described 
by the Schwarzschild metric, and we can start study- 
ing them by investigating the properties of this met- 
ric. One of the most important things that one should 
study first is the classical stability of a given solution. 
In fact, if a solution is not stable, then it will most 
certainly not be found in nature, unless the instabil- 
ity timescale is much larger than the age of our uni- 
verse. What does one mean by stability? In this classical 
context, stability means that a given initially bounded 
perturbation of the spacetime remains bounded for all 
times. For example, the Schwarzschild spacetime is sta- 
ble against all kinds of perturbations, massive or mass- 
less ij. Thus, the Schwarzschild geometry is indeed ap- 
propriate to study astrophysical objects. On the other 
hand the Kerr spacetime, describing a rotating black 
hole, is stable against massless field perturbations but 
not against massive bosonic fields (although the in- 
stability timescale is typically much larger than the age 
of the Universe, thus presenting no real danger). In four 
dimensional (asymptotically flat) spacetime, the Kerr- 
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Newman family is the most general black hole solution 
to Einstein equations, and, if we exclude massive bosons, 
they are all stable. However, unstable solutions seem to 
be more common than previously thought. Take for ex- 
ample anti-de Sitter (AdS) spacetimes. This has become 
a very popular background spacetime since it was conjec- 
tured [Sj that there is a duality between the gravitational 
degrees of freedom in the bulk of AdS space and a Con- 
formal Field Theory formulated on the boundary of that 
space; this is the AdS/CFT correspondence conjecture. 
A black hole in this spacetime corresponds to a thermal 
state on the GET. Are AdS black holes stable? Not all, 
as shown recently by Gardoso and Dias 3, who have 
proved that small Kerr-AdS black holes are classically 
unstable. This instability is due to a "black hole bomb" 
effect '3], whereby waves are successively amplified near 
the black hole event horizon and reflected at the bound- 
ary of the AdS spacetime (the special thing about AdS 
spacetime is that its boundaries, spatial infinity, works 
as a wall) . If we now consider higher dimensional space- 
times, which are of interest to string theory and/or extra 
dimensional scenarios, instabilities seem to be much more 
common. For instance, even though higher dimensional 
Schwarzschild black holes 0,0 are stable @, their rotat- 
ing counterparts seem not to be, at least for large rota- 
tion. Indeed, it was proved by Gregory and Laflamme 
that black branes are classically unstable against a sec- 
tor of gravitational perturbations (the tensorial sector) , 
and this result was used recently by Emparan and My- 
ers jlO| to argue that ultra-spinning higher dimensional 
black holes should be similarly unstable (recall [jj that 
for spacetime dimensions D greater than 5, -D > 5, there 
is no limit on the rotation parameter). Recently Gardoso 
and Lemos [TH have uncovered a new universal instability 
for rotating black branes and strings, which holds for any 
massless field perturbation. The gist of their argument 
is that transverse dimensions in a black brane geometry 
act as an effective mass for the fields, which simulates a 
mirror enclosing a rotating black hole, thereby creating 
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a black hole bomb 0| . 

For other types of instabilities see for example 
Here, we shall investigate the stability against scalar 
perturbations of ultra-spinning black holes. In four and 
five dimensions, there is an upper bound for the rotation 
parameter a of a Kerr black hole , and when the black 
hole saturates that bound we say it is an extremal black 
hole. Now, it is known J^sil^ JfiJ that the characteristic 
frequencies (quasinormal frequencies, or QN frequencies) 
of four or five dimensional Kerr black holes always have a 
negative imaginary part (the field is decomposed accord- 
ing to ^' ~ e~"^*$(r, angles)) so the Kerr spacetime is 
stable. However, as the black hole approaches extremal- 
ity, the imaginary part of the QN frequency tends to 
zero, thus raising the possibility that if there was no up- 
per bound on a the QN frequencies could have a positive 
imaginary part, or in other words, the spacetime could 
be unstable. This will be our main motivation for this 
study. We are not studying tensorial perturbations, so 
we shall not be dealing with Gregory-Laflamme type of 
instabilities. Instead, we are more interested in finding 
out what are the consequences, if any, of having arbitrar- 
ily large angular momentum for a black hole. We shall 
focus, for concreteness, on six-dimensional rotating black 
holes, but we suspect that the general features born out of 
this study are valid for any spacetime dimension greater 
than five. Previous work on related subject includes that 
of Ida et al and Berti et al 0| who studied five 
dimensional Kerr black holes, in large and compact ex- 
tra dimensions respectively. As we remarked, in five di- 
mensions there is a bound on the rotation of the black 
hole, and thus these works are not suitable for studying 
possibly new phenomena appearing for unbound rotation 
parameter. 
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and dfl^ denotes the standard metric of the unit n- 
sphere. This metric describes a rotating black hole in 
asymptotically fiat, vacuum space-time with mass and 
angular momentum proportional to /i and fj,a, respec- 
tively. Hereafter, /z, a > are assumed. 

The event horizon is located at r = rn, such that 
A|r=rff = 0, which is homeomorphic to 5'^+". For n = 0, 
the standard 4-dimensional case, an event horizon exists 
only for a < /i/2. When n ~ 1, an event horizon exists 
only when a < ^JJi, and the event horizon shrinks to zero- 
area in the extreme limit a — > ^JJi. On the other hand, 
when n > 2, which is the part of the parameter space 
which we shall focus on, A = has exactly one positive 
root for arbitrary a > 0. This means there is no bound 
on a, and thus there are no extreme Kerr black holes in 
higher dimensions. 



B. Separation of variables and boundary conditions 

Consider now the evolution of a massless scalar field 
^I^ in the background described by Q). The evolution is 
governed by the curved space Klein-Gordon equation 



0, 



(4) 



II. FORMULATION OF THE PROBLEM AND 
BASIC EQUATIONS 

A. The background metric 

Here we adopt the notation of Ida et al 16] , and we also 
correct some typos appearing in their equations. In four 
dimensions, there is only one possible rotation axis for 
a cylindrically symmetric spacetime, and there is there- 
fore only one angular momentum parameter. In higher 
dimensions there are several choices of rotation axis and 
there is a multitude of angular momentum parameters, 
each referring to a particular rotation axis Q- Here we 
shall concentrate on the simplest case, for which there is 
only one angular momentum parameter, which we shall 
denote by a. The metric of a (4 -|- n)-dimensional Kerr 
black hole with only one non-zero angular momentum 
parameter is given in Boyer-Lindquist-type coordinates 

by 



A - sin^ -d 



-dr 



2a(r2 + - A) sin^ i? 



dtd(f 



where g is the determinant of the metric. The metric 
appearing in Q should describe the geometry referring 
to both the black hole and the scalar field, but if we con- 
sider that the amplitude of ^E* is so small that its contri- 
bution to the energy content can be neglected, than the 
Kerr metric (^Q should be a good approximation to g^^^, 
in Q . We shall thus work in this perturbative approach. 
It turns out that it is possible to simplify considerably 
equation Q if we separate the angular variables from 
the radial and time variables, as is done in four dimen- 
sions This separation was accomplished, for higher 
dimensions, in |l9l | for five dimensional Kerr holes (who 
work with two spin parameters) and also in (2^ for a gen- 
eral 4+n-dimensional Kerr hole. Since we are considering 
only one angular momentum parameter, the separation 
is somewhat simplified, and we can follow |16[. In the 
end our results agree with the results in |19l l2Cll | , if we 
consider only one angular momentum parameter in their 
equations. 

We consider the ansatz cj) = e''^^'''^^ R{r)S{d)Y {Vl), 
and substitute this form in Q), where Y{Vl) are hyper- 
spherical harmonics on the n-sphere, with eigenvalues 
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given by + n — 1) (j = 0, 1, 2, • • •). Then we ob- 

tain the separated equations 



1 f d ■ a n ad's 

- — — sln^9cos"^?— 

Sl^ i} cos" ■& yd-d d-d 



ij?a^ cos^ -d 
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z9-j(i + ?^ -l)sec2i9-|-^] 5 = 0, (5) 



and 



dr \ dr 



[cj(r^ + a^) — ma] 



-\'>R = 0, 



(6) 



where X := A ~ 2'muja + uP'c?' . 

The equations (O and 10 must be supplemented by 
appropriate boundary conditions, which are given by 



(r - thT 



as r ^ rn , 
as r ^ oo . 



(7) 



where 



[{r'jj + a^)ijj — too] th 



(8) 



has been determined by the asymptotic behavior of the 
Eq. jnj. In other words, the waves must be purely ingoing 
at the horizon and purely outgoing at the infinity. For as- 
signed values of the rotational parameter a and of the an- 
gular index I ,j ,m there is a discrete (and infinite) set of 
frequencies called quasinormal frequencies, QN frequen- 
cies or loqn, satisfying the wave equation © with the 
boundary conditions just specified by Eq. Q. The QN 
frequencies are in general complex numbers, the imagi- 
nary part describing the decay or growth of the perturba- 
tion, because the time dependence is given by e~'"*. We 
expect the black hole to be stable against small pertur- 
bations, and therefore ujqn is expected to have a positive 
imaginary part, so that the perturbation decays exponen- 
tially as time goes by (recall that the time dependence 
of the wavefunction is e*'^'). As usual, we will order the 
QN frequencies ujqn according to the absolute value of 
their imaginary part: the fundamental mode (labeled by 
an integer n = 0) will have the smallest imaginary part 
(in modulus), and so on. We refer the reader to |2j and 
|22i] for further details on QN frequencies. 



III. NUMERICAL COMPUTATION 

Now that we have a well posed problem, we have to 
solve for the characteristic QN frequencies. The most 
powerful method to date is that of Leaver jQ, which 
makes use of a continued fraction representation, and 
which can determine the resonant frequency lu and the 
separation constant A with very high accuracy. We as- 
sume the following series expansion for S 



S = (sin7?)l"l(cosz9)J'^afe(cos2z9)'^ 



(9) 



which automatically satisfies the regular boundary con- 
ditions at i} = 0,7r/2 whenever converges. Substituting 
this into Eq. 0, we obtain the three-term recurrence 
relations 

aofli + Poao = , 

akak+i + hak + IkO-k-i = , (A: = 1, 2, • • •)(10) 

where 

Uk = -2(A:+l)(2j + n + 2/2 + 1), 

Pk = {] + \m\ + 2k){j + n+\m\ + 2k + l) - A, 

2 2 

Ik = -LU^a^ . 

Here, we have defined the dimcnsionless quantity := 
Lurn and a* :— a/rn, since the behavior of the system 
depends only on a*. When a* = 0, the eigenvalue A 
is explicitly determined from the requirement that the 
series expansion ends within finite terms, since otherwise 
divergent. Thus we have 

A = {2£ + j+\m\){2e + j + \m\+n + l) + 0{uj^a*), 
{£ = 0,1,2, ■■■) (11) 

and the Oth-order eigenfunctions are given in terms of 
the Jacobi polynomials: 



P,,„ = (sint?)l™l(cos?9)^' 

xF(~£,£ + j + \m\ + 



n+ 1 , n+ 1 ^ 
— : cos -u 



(12) 



In a similar way, for the n > 1 case, we expand the 
radial function R into the form 



R 



r-rn 



OC 



r 



k 



-(n+2)/2-i(T 



(13) 



where 6o is taken to be fog = 1- If n = 2, the expansion 
coefficients hk in equation (|13|l are determined via the 
seven-term recurrence relation (it's just a matter of sub- 
stituting expression H13|l in the wave equation ©), given 
by 

aoh + f3obo = , 
aib2 + Pibi -1-7160 = 0, 
5263 + P2b2 + 7261 + ^260 = , 
5364 + Psbs + 7362 + Ssbi + £360 = , (14) 
5465 -I- P4b4 + 7463 + (^462 -I- €461 + C4^o = , 
akbk+i + Pkbk + Ikbk-i + Skbk~2 + hbk-a 
+Ckbk-4: + Vkhk-b ^ Q , (fc = 5, 6,---) 



where 



fc=o 



Ctk 



(l + fc)(l + fc + 2^^7)(3 + <)^ 
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/3fe 



-18 
9fc(3 



Ik 



efe 



Cfe 



36fc2 
- 8i(T H 



3A - 27ia + 36(7^ + IScrw, + 2wJ 



{12 + 3j(i + 1) + 30fc2 + 2m2 + A + 21icr 
-30ct^ + 3fc(7 + 20icr + Aiu^) + Qiuj^ - 12ctw4 
at {2 + j{j + 1) + 6fc2 + 4icr - 6(7^ 
+2fc (2 + 6i(T + icj,) + icj, — 2(Tijj» + 2uj1} , 
-2mLu^al + at{7 + + 1) + ISfc^ - lOicr 

2icj, — Bcrti)* + 2ik{5i + 15a + Alu^)} 
A — llicr — 20(7^ — Siujtf 



3{9 + 20fc2 

-2tj? 



ifc(lli 



-40fT 

an30 + 9j(i + l) + 62fc2 



12tj,)} 
+ 2A 



- 12aa;* 
38ia 



—62a^ — 9iujr — Sdauj^ — Guil 

+ 2ifc(19i + 62cr+ 18w,)} , 
= -66 - 54fc2 - A + lQ2ia + 54(7^ 
+ 6fc(17 — 18i(7 — 5iw») + 30iw» + SOauj^ + AluI 

- 2at {14 + 3j{j + 1) + lOfc^ - 20icr - lOcr^ 

— 6i[x), — Gcto;, — ijjJ + 2ik{10i + 10a + 3ijJ*)} 

- al {90 + 10j{j + 1) + 68fc2 + A - 132icr 

-68(7^ - 40itJ, - 40ctw, - 6ljI 
+4ik{33i + 34cr + lOcj,)} , 
28fc^ - 91i(7 - 28(7^ - 21iw» - 12auj^ 
12uj.) 



81 

+ ifc(9ii + 56(7 
a*{47 + 4j(j + l) 

— 14zcj* — Scrcj^ — - 



15fc" 



4w,)} 



50icr - 15a^ 
2ik{2i + 15(7 4 
al {126 + 5jO' + 1) + 42fc2 - 138i(7 - 42a2 
—35iuJt: — 20(7w, — 2lo\ 
+ 2ifc(69i + 42(7+ lOw,)}, 
-(1 + a^)[44 + 8fc2 - 37i(7 - 8ct^ - 5iw, - 2(7w 
+iA;(37i + 16(7 + 2a),) 

{34 + j(j + 1) + 6fc2 - 28i(7 - 6(7^ 
— 5ia;* — 2(7ijj* + 2ifc(14i + 6ct + cj*)}] , 
(l + a2)2(-3 + fc + i(7)2. 



By making a Gaussian elimination four times, one can 
reduce the seven-term recurrence relations (I15|l to the 
three-term recurrence relations, which is given by 



, A: = l,2,.. 



(15) 



For more details on how to obtain the coefhcients of the 
three-term recurrence relation, we refer the reader to |23j | . 
An eigenfunction satisfying the QN mode (QNM) bound- 
ary conditions behaves at the event horizon and infinity 
as Eq. 0. Therefore, one can see that the expanded 
wave function (|13() satisfies the QNM boundary condi- 
tions if the expansion in (|13|) converges at spatial infin- 
ity. This convergence condition for the expansion (|13|l . 
namely the QNM conditions, can be written in terms of 



TABLE I: In this Table we compare our numerical results 
for the QN frequencies of six-dimensional, non-rotating black 
holes, with results obtained through WKB techniques, and 
we also indicate the error involved using the WKB approach. 
The results refer to the fundamental mode of several l,j,m 
perturbations. 
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, ,Num„ 
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%Im 











0.8894-f0.5331i 


0.7682-f0.5265i 


13.6 


1.2 
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1.4465-f0.5093i 


1.3846+0.4933i 


4.3 


3.1 


1 





1 


2.57914-0.49891 


2.5455+0.4942i 


1.3 


0.9 


1 


1 


1 


3.14784-0.49731 


3.1205-f0.4944i 


0.9 


0.6 



the continued fraction as [13. l2^ 



/3o 



a'oli «i72 5273 

K- 02- h- ' 



ll2 



0,(16) 



/32- 



where the first equality is a notational definition com- 
monly used in the literature for infinite continued frac- 
tions. Here, we shall adopt such a convention. As for 
the determination of the separation constant A, exactly 
the same technique of the continued fraction can be ap- 
plied. The continued fraction equation for the separation 
constant is then given by 



q;o7i ai72 "273 _ „ 



(17) 



In order to obtain the QNMs, one has to solve numeri- 
cally the two coupled algebraic equations (|16|l and H17|l . 
following for example the procedure in pSil] . 



IV. NUMERICAL RESULTS 

Using the method described above, we have made 
an extensive search for the QN frequencies of six- 
dimensional rotating black holes, for several values of 
the parameters I , j ,m. The numerical results are sum- 
marized in Table HI and in Figs. I1I2I To check our 
code we have first computed the QN frequencies of six- 
dimensional Schwarzschild black holes, and compared 
them with analytic WKB results [2^. The results are 
shown in Table ^ along with a computation of the er- 
ror involved in the WKB calculation. Although for 
I = j = m ^ the error is quite large (14% for the 
real part) it quickly decreases as I ,j ,m, increase. We 
can thus say that the code is tested (or then that the 
WKB approximation yields good results...). The full nu- 
merical results for six-dimensional rotating black holes 
is presented in Figs. I1I2I In Fig. we present the re- 
sults referring to the real part of the fundamental QN 
frequency, as a function of a, for several I ,j ,m values. 
The real part ofuQN seems to decrease monotonically as 
a/rn increases, and for very large a/rn, it asymptotes 
to zero. 
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FIG. 1: Real part of the fundamental QN frequency as a 
function of the rotation parameter a for some I ,j ,m values. 
The maximum is reached at zero rotation, and as a increases 
the real part of wqjv decreases monotonically. 
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FIG. 2: Imaginary part of the fundamental QN frequency as 
a function of the rotation parameter a for some I ,j ,m values. 
Notice that, for all values of a the imaginary part is always 
positive, -which means that even ultra-spinning black holes are 
stable. 



In Fig. |21we present the results referring to the imagi- 
nary part of the fundamental QN frequency, as a function 
of a, for several I , j ,m values. Although the pattern is 
more complex now, one can see that up to a/rn — 40 
the imaginary part of loqn is still positive, and thus the 
modes are stable. It proves very difficult to get numerical 
results for higher values of a/rn, but if the trend contin- 
ues, and we have no reason to believe otherwise, it looks 
like an instability will never set in. Results for higher 
overtones, both for the real and imaginary part follow a 
similar trend. 



V. ANALYTICAL RESULTS 

We concentrate on n > 1 (six and higher dimensions). 
Let us first re-write the radial wave equation in terms of 
dimensionless variables. 



y-"(y"Ai?')'+[ 



A 



y 



A + 2mLu^a^-Luiai\R = 0, (18) 



where y = r/rn and A = A/rj^ — + aj — (1 + aj)j;^ 
The horizon is at y = 1 which is the real root of A = 0. 

In order to bring the radial wave equation into a 
Schrodinger-like form, it is convenient to introduce 



= y"/^(5(i/))'/'i?(y) , 9{y) - iv' + a^A 

In terms of 5", the wave equation (|18|l reads 

- h{y) {h{y) *')' + V[y)-^ = [u, - mn.f-^ 
where 

A 



(19) 
(20) 



Ky) 



V9¥) 



(21) 



The potential is 

9{y) 



V{y) = A- 2mLu^n^ - m^ajrj 



2„2n2 , rL{n + 2) 



+ n-l) + n/2(n/2 - l)]a'^ n'^{l + a^) 



Ay 



n+l 



A 



[(mf7, — 2a;»)(y^ + a^) -f ma] + 
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9 



and we have introduced the angular velocity of the hori- 
zon, 



H 



(23) 
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The potential vanishes at the horizon {y = 1) and ap- 
proaches a constant at infinity {V — > mf2»(mr2* — 2w*) 
as 2/ — > oo). In general, it depends on the frequency oj. 
However, in the two extremal limits a — *■ (Schwarzschild 
limit) and a — > c» (where one expects an instability to 
develop), the dependence on lo drops out due to the van- 
ishing of the angular velocity [Q, 0) in these limits. 
In terms of the tortoise coordinate defined by 



dy 
dy* 



Ky); 



the wave equation H20|) may be written as 



^ 2 + y\v{v*)] * = (w* - mn^f^ . (25) 

dyi 

Let us consider the two extremal limits separately. In the 
limit a — > 0, we obtain the Schwarzschild wave equation 



dyl 



where 



(n + 2)2 



Ay 



ri+3 



(26) 



(27) 



and L = 21+ j + \m\ + ^-!-^ . To estimate the eigenfrequen- 
cies, expand around the maximum, ymax of the potential. 
Setting V^iy 

max 

) = 0, we obtain 



n + 3 



l/(n+l) 



+ o{l/L). (28) 



Expanding around the maximum, we may approximate 
the potential by 



Voiyiy.)] « - p^y, ~ y.{y„,ax)? , (29) 



where 

and 

13' = 



(30) 



2 dyl 



^i^iymax)) ^0 {ymax) ■ 



2/.=a.(amax) 



Explicitly, 



n -I- 1 



n -I- 3 V 'T- + 3 



2/(n+l) 



L' + Oil); 



(31) 



(32) 



0' 



(n + l)= 



4/(n+l) 



L^+o{l). (33) 



(n + 3)2 V" + 3, 
The wave equation becomes 
_^"_/32^2^^ (^2_^2)^^ x = y,-y,(2/™,,). (34) 



The solutions obeying the right boundary conditions at 
X —^ ±oo are 



^Hj^ = HN{\fiPx)e'^''^/^ , iV = 0,l,2,, 



(35) 



where are Hermite polynomials, with corresponding 
eigenvalues 



(36) 



(24) Exphcitly, 



cj, = C(n) {i + i\/^^TT(iV+ i)} +o(l/L), (37) 



with 



C{n) 



(38) 



This result is exactly what one gets by using a standard 
WKB approach 

Turning to the other extreme, a — > oo, we have 



h{y) w y 



1 - 

giv) 

(«-l)/2 _ . 



1 



1 



(n-l)/2 



y 



(39) 
(40) 
(41) 



and the potential becomes to leading order in 1/a, 

1 



J + 



n-1 



Vooiy) 

2 



= y 

n — 3 



,n — 3 



1 



y 



(n + 1)2 

IQyn-l 



(42) 



where we assumed A < o(a^). The wave equation has a 
well-defined limit as a — > oo. However, the potential is 
positive and diverges as y — + cxd for n > 3, so sublead- 
ing terms are needed to estimate the eigenfrequencies. 
For n < 3 (six and seven dimensions), u; approaches a 
constant value independent of a which is easily found by 
solving the Schrodinger equation. This asymptotic value 
only depends on j. 

In six dimensions (n = 2), the potential exhibits a 
maximum and may be approximated by an inverted har- 
monic oscillator potential, as in the Schwarzschild limit. 
The frequencies can be found explicitly as functions of 
j taking advantage of the fact that the equation for the 
maximum V^{ymax) = is quadratic. In Table ITTl we list 
the QN frequencies as a function of j, obtained using this 
analytical scheme for a oo. 



VI. CONCLUSIONS 

We have investigated numerically the stability of six- 
dimensional rotating Kerr black holes, with one rotation 
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TABLE II: In this Table we show the results of an analytical 
WKB type scheme for computing the QN frequencies in the 
ultra-spinning regime, a oo. The results depend only on j. 
This scheme shows that uqn asymptotes to a constant value, 
which is consistent both qualitatively and quantitatively with 
the numerical results shown in Figures and |5| 



J 



Analy 



+ 0.1621 



0.576 + Oi 



1.078 + Oi 



they suggest that the Gregory-Laflamme [9| instabihty 
should be the cause). This, for the moment, is a major 
challenge specially because there is no known formalism 
to handle gravitational perturbations of higher dimen- 
sional Kerr black holes. Such a formalism could also 
prove useful in studying at depth the recently discovered 
instability for rotating black branes and strings 
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